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In this paper the charged Higgs signal through the decay to a pair of muon and neutrino {H^ — >■ 
^u) is analyzed. The analysis attempts to estimate the amount of muonic signal of the charged Higgs 
at LHC at a center of mass energy of 14 TeV. The signal process is the top quark pair production 
with one of the top quarks decaying to a charged Higgs (non SM anomalous top decay) and the 
other decaying to a W boson which is assumed to decay hadronically to two light jets. Due to the 
small branching ratio of charged Higgs decay to muon, results are quoted for data corresponding to 
an integrated luminosity of 300 fb^^ which is expected to be collected at the LHC high luminosity 
regime. It is shown that a signal significance close to 5(t down to below Icr is achievable for a charged 
Higgs mass in the range 80 GeV < m{H^) < 150 GeV taking the top quark pair production with 
both top quarks decaying to W bosons as the main irreducible background. 

PACS numbers: 



I. INTRODUCTION 

The Large Hadron Collider (LHC) at CERN is one of the biggest high energy physics experiments which is expected 
to provide fruitful information about the possible new phenomena at the scale of sub-atomic collision events. Currently 
running at a center of mass energy of 7 TeV, LHC has delivered almost 5/6"^ data to each one of the CMS and 
ATLAS experiments. Obviously one of the exciting searches is the search for the Higgs boson which is expected to be 
responsible for the mass of particles through the Higgs mechanism [l-^Sj . The current results of Higgs boson searches 
in CMS exclude the Higgs boson mass in the range 127-600 GeV at 95% confidence level [1], while the excluded ranges 
by ATLAS are 112.7-115.5 GeV, 131-237 GeV and 251-468 GeV at 95% confidence level 0. Both above resuhs 
obviously confirm and extend the current Tevatron results [8] . 

The above news imply that the next few years of LHC experiments will be very exciting and provide a variety of 
information which develop our knowledge and level of understanding of the Standard Model and in case of observing 
the Higgs boson a reasonable understanding of the origin of mass spectrum through the Higgs mechanism is achieved. 
While an extensive effort is ongoing at the Tevatron and LHC experiments in search for the Standard Model Higgs 
boson, there are arguments that a single Higgs boson predicted through the SM Higgs mechanism is not a complete 
theory. In fact the problem of quadratic divergence of the Higgs boson mass when quantum corrections are included 
does not receive a solution as elegant as the Supersymmetry (SUSY) which introduces a super-partner for each particle 
whose spin is different by half unit ^2] . Making a double universe this way, each fermion acquires a boson as a partner 
and vice versa. It is shown that within these models the quadratic divergence of the Higgs boson mass reduces to a 
logarithmic function of the cut-off of the theory. 

One of the interesting results of supersymmetric models is that instead of a single Higgs boson, several such particles 
with different properties are produced. In the minimal supersymmetric extension to Standard Model, the so called 
MSSM, five physical Higgs bosons are predicted. Among these five, two are charged and the other three are neutral. 
It might look difficult at the beginning to distinguish between the SM neutral Higgs boson and a neutral MSSM Higgs 
boson. However the decay preferences of the SM Higgs boson and MSSM neutral Higgs bosons are different, i.e. there 
are different branching ratio of decays when SM Higgs and MSSM neutral Higgs bosons are compared and therefore 
different strategies are taken in the search for them [l^, . 

A long time attempt has been dedicated to the search for the MSSM neutral Higgs bosons in the past and current 
experiments. The LEP results on such searches are reported in [l^ where the lower limits on the mass of neutral 
Higgs bosons are found to be m^jo > 91.0 GeV and^m^o > 91.9 GeV and low tan/3 values, i.e. 0.5 <tan/3 < 2.4 are 
excluded. The Tevatron results are also quoted in [l4j where high tan/3 values are excluded. The LHC collaborations 
have also started an intensive search for such neutral MSSM Higgs bosons. The CMS Collaboration has recently 
published the first results of such searches in [l^ which has set the strongest limits on tan/3. 

Another direction of search is looking for the charged Higgs. Being a charged particle, it provides unique signatures 
which makes it distinguishable from neutral Higgs bosons. The decay channels are different and due to being a 
spinless particle, it produces different topology and kinematics of events when compared to the electroweak gauge 
bosons such as . Observation of such a different particle is a crutial signature of models beyond Standard Model 
such as MSSM. This is why a special care has been taken in searches for this particle. 
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The LEP search for the charged Higgs has been reported in [To!] where a charged Higgs mass below 78.6 GeV is 
excluded at 95% confidence level. This limit is of course softer than that obtained in indirect searches in [l^ where a 
charged Higgs mass below 125 GeV is excluded at high tan/3 values (tan/3 > 10). The Tevatron searches are reported 
in [la-il by the DO Cohaboration and [H-i^ by CDF CoUaboration. The overall result is 2 < tan^ < 30 for 
m{H^) > 80 GeV while wider tan/3 values are available for higher charged Higgs masses. It should be noted that 
different strategies are used in the search for this particle. In addition to direct search limits and limits obtained from 
searches for the neutral MSSM Higgs bosons, indirect search limits are obtained using analyses which look for any 
deviation in the measured SM cross sections as a hint for the existence of beyond SM particles such as the charged 
Higgs bosons. A review of direct and indirect searches at the Tevatron can be found in [28| . 

Results from B-factories are also used to put some constraints on the charged Higgs mass as the existence of a charged 
Higgs in B-meson decay diagrams may alter branching ratio of decays and can be verified experimentally. A study of 
the & — > S7 transition process using CLEO data has excluded a charged Higgs mass below 295 GeV at 95 % C.L. in 
2HDM Type II with tan/3 higher than 2 in [2^. This limit is obviously stronger than direct search limits, however, 
there have been studies such as the one reported in [sO] where it has been shown that in a general 2HDM with complex 
Higgs-fermion coupling, the charged Higgs mass can be as low as 100 GeV. The above direct searches carried out 
at LEP, Tevatron and LHC are based on the charged Higgs decay to tau and the tau neutrino, i.e. 7?+ — T~^i/r- 
This is because the light charged Higgs decays predominantly to the tau lepton and its signature is the excess of tau 
leptons over the Standard Model expectation which is also inferred as some imbalance in the number of leptons from 
W decays. This is the so called lepton universality breaking phenomenon. The charged Higgs is produced from an 
on-shell (off-shell) top quark if it is lighter (heavier) than the top quark. 

The r lepton identification algorithm follows theoretical and phenomenological studies of different decay channels of 
this particle. These studies have improved the r lepton identification over the years. An example of a comprehensive 
theoretical study of the t lepton decay channels can be found in [l^. In [3| and [13] algorithms for the r lepton 
identification were proposed for Tevatron and LHC. In a different work an update of the algorithm specific to the case 
of heavy charged Higgs was proposed in [l^ . The charged Higgs analyses searching for the tv decay mode extensively 
use the idea proposed in above references. 

The main process for the light charged Higgs is tt production with at least one of the top quarks decaying to a charged 
Higgs, i.e. tt -> H^W^bb. This is to be compared with a background from tt — ^ W^W^bb. Therefore the top decay 
to charged Higgs, i.e. t — >■ H^b is compared with t — W^b and since W bosons equally decay to electron, muons and 
taus, the existence of a charged Higgs in the top decay is concluded by observing an excess of tau leptons over what 
we expect from t W'^b. 

A helping character of top decay to charged Higgs is that the mass of the charged Higgs is taken to be higher than the 
W boson mass required by the direct search results of LEP reported in [l^ . Therefore a more energetic tau lepton 
produced in charged Higgs decay is better identified than its partner from W boson decay. The same arguments apply 
to a heavy charged Higgs which is produced in gluon-gluon interactions producin g a top quark in association with 
the charged Higgs, i.e. gg — >■ tbH~ . The details of this process can be found in (sTI. |32||. Again for the heavy charged 
Higgs, the main background is taken to be the Standard Model process tt — W^W~bb and the same approach as 
above is used for the search. 

In CMS and ATLAS collaborations, several Monte Carlo studies have been carried out for the light and heavy charged 
Higgs. The CMS studies of the light and heavy charged Higgs decay to tau are documented in [H, 113]. The heavy 
charged Higgs decay to top and bottom quarks has also been studied in CMS i n fSSll. Although the latter does not 
lead to a significant signal due to the large hadronic background, the analyses in [33,[3J] clearly show that a large area 
of parameter space in terms of m{H^) and tan/3 is available for observation of this particle. The ATLAS collaboration 
also quotes their search results on the potential of observing a charged Higgs boson in different final states in [l^ . 
An overall conclusion by CMS and ATLAS collaborations is that the charged Higgs decay to tv will be the most 
promising decay channel in the search for this particle and will certainly be taken as the highest priority channel for 
the early data analysis at LHC. That is why no more attention has been paid to the charged Higgs rare decays such 
as decay to muon in the current searches. Since — >■ tv provides almost enough statistics for a wide range of the 
charged Higgs mass, any possible discovery is expected to happen in — >■ tv channel. 

In this paper the attempt is to verify the observability of charged Higgs decay to muon quantitatively. The aim is 
to verify whether this signal would be finally observable at the LHC high luminosity regime or not. To the best of 
our knowledge, no detailed analysis has been published in the literature focusing on this channel. Obviously this 
channel is not suitable as an early data analysis and can not be considered as a discovery channel because if a charged 
Higgs exists, the tauonic decay signal would have much more statistics to tell us about the existence of the charged 
Higgs particle. However by the time LHC will have provided enough data to study the rare decays of this particle, 
the muonic signal would be a good candidate providing a clean signature for studies of the charged Higgs boson 
spin, coupling to fermions, etc. For instance, studies of the muon fiight direction and pt distribution could provide 
opportunity to verify the charged Higgs spin through the shape fits in a similar way as has been done for bosons 



3 





Signal 


Background 




80 GeV 


100 GeV 


120 GeV 


150 GeV 




o-xBR 


302 fb 


228 fb 


149 fb 


38 fb 


124 pb 



TABLE I: Signal and background cross sections times branching ratios. For the signal tan/3=20 has been used. 



[36| . A detailed study of such phenomena could be carried out by LHC experiments when enough data is available. 
As a comparison, the charged Higgs decay to muon suffers from a small coupling which is proportional to the lepton 
mass and the branching ratio in this case is roughly two orders of magnitude smaller than the case of tauonic decay. 
This is compensated to a reasonable extent by the fact that muons are identified with a high efficiency thanks to the 
good detector functionality and the powerful muon detection systems of LHC experiments. On the contrary, the tau 
lepton identification requires a sequence of selection cuts and the total efficiency is low [37] . Therefore the charged 
Higgs signal in this case is reduced significantly from what is produced at the beginning although still being higher 
than that of charged Higgs decay to muon. The decay to c6, cs or us does not lead to a promising signal as they all 
suffer from the large hadronic background from QCD multi-jet events. 

Above thoughts motives a study of the charged Higgs decay to muon in a quantitative way and estimation of the 
signal significance as a function of the charged Higgs mass. Since the light charged Higgs comes mainly from the tt 
production which has a high cross section at LHC, this study focuses on the light charged Higgs which is expected 
to give the highest sensitivity. Extrapolation to the case of heavy charged Higgs can be made later or in a separated 
study. Therefore this is a study with the goal of evaluating the possibility of observing in the region of 

light charged Higgs, i.e. below the top quark threshold. 



II. SIGNAL AND BACKGROUND IDENTIFICATION AND SIMULATION 

The signal for this analysis is tt — H^W^bb followed by ^iv and jj where only the hadronic decay 

of the W boson is considered to be able to reconstruct the Higgs boson transverse mass. Therefore the final state 
consists of a muon, four jets two of which are b-jets, and some amount of missing transverse energy originating from 
the neutrino in the event. Figure [T] shows the top quark branching ratio of decay to W and charged Higgs as a function 
of the charged Higgs mass. The charged Higgs masses under study in this work start from the direct search lower 
limit in [l^ and cover a mass range up to 150 GeV where the signal cross section becomes very small. The tan/3 used 
in the study is 20, as higher values are close to or inside the excluded area in [l^l- The charged Higgs branching ratio 
of decays are shown in Fig. [2] As is observed the decay to muon has a branching ratio of roughly 0.35 x 10~^. The 
lighter charged Higgs, the higher branching ratio of top decay to charged Higgs. However it is below 0.1 for all the 
region of study and when multiplied by the charged Higgs branching ratio of decay to muon results in a small cross 
section compared to background. 

The main background is ti — W^W^bb followed by — > /ii/ and — jj. Table |T] lists the signal and background 
cross sections multiplied by branching ratios. In calculating the numbers quoted in Tab. U MCFM 5.7 has been 
used with the parton distribution function (PDF) set to MRST 2006 which was taken by linking LHAPDF 5.8.3 
(39| to MCFM. The top quark mass is set to 175 GeV. The total ti cross section using the above PDF set is 879 
pb. The W boson branching ratios are taken from the Particle Data Group (PDG) data [ioj (BR(T4^ — >■ jj^ 



0.68 and BR(Vt^ — iiv) ~ 0.106). The charged Higgs branching ratios are calculated using HDECAY 3.51 41 1. 
The event simulation and calculation of branching ratios of top quark decays are done using PYTHIA 6.4.21 42 1. 
When generating events the same PDF set is linked to PYTHIA through the LHAPDF interface. The study is 
based on a parton level generation of events followed by hadronization and jet reconstruction using the PYTHIA jet 
reconstruction tool PYCELL. The initial and final state QCD and QED radiations are also taken into account in the 
simulation settings. The jet reconstruction cone opening angle is 0.5 and only jets with pseudorapidity \r]\ < 3 are 
reconstructed. Here the pseudorapidity is defined as rj = —ln{ta.n{9/2)) where 6 is the polar angle with respect to 
the beam axis. The MSSM, TO/j-max scenario is used with the following parameters: M2 = 200 GeV, Mg = 800 GeV, 
At = 200 GeV and Msusy = 1 TeV. 



III. EVENT ANALYSIS AND SELECTION 



In this section signal and background events are generated and their kinematic properties are compared and a set 
of selection cuts are introduced for the signal selection. Since the signal cross section is small, the analysis is designed 
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FIG. 1: Branching ratio of top decay to charged Higgs and W boson as a function of its mass. 
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FIG. 2: Branching ratio of charged Higgs decays a function of its mass. 

for high himinosity regime of LHC. Therefore a total integrated luminosity of 300 fb^^ is used to obtain the first 
estimate of the signal size. 

Since signal events produce a muon, two light jets, two b-jets and some missing transverse energy through the 
following chain tt — >■ H^W^bb — >■ fii'jjbb, events can be triggered online by requiring a single muon trigger. Figure[3] 
shows the muon transverse momentum distributions in signal events with varying charged Higgs mass as well as the 
background distribution from the SM top pair production where the muon comes from the top decay to a W boson 
through tt W^W^bb nvjjbb. 

This is a mass independent analysis, i.e. the signal selection cuts do not depend on the charged Higgs mass. A 
more detailed study by variating selection cuts as a funcation of the charged Higgs mass may produce higher signal 
significances. However that requires a different study including cuts optimization which can be performed in a 
separated study. In the following a common set of selection cuts are applied on all charged Higgs mass cases. The 
cut on the muon pT and |?7| is applied as the following: 
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FIG. 3: Muon transverse momentum distributions in signal events with difTerent charged Higgs massed and the background 
sample. 



Muon pT > 30 GeV, \tj\ < 2.5 (1) 

The above requirement selects muons in the central detector region with the momentum threshold above 30 GeV. 
Harder cuts are avoided due to the small cross section of the signal, and the fact that they do not improve the overall 
significance of the signal. This was verified by evaluating t^gj where £5(6^) is the signal (background) relative 
efficiency when the \th cut is applied. This ratio is a measure of goodness of a cut when the signal significance is 
going to be calculated as S/^/B. The next step starts with counting number of jets in the central region of \rj\ < 2.5 
and the jet transverse energy higher than 30 GeV. Softer jets (with lower Et) are not taken into account to avoid 
reconstruction inefficiency and jet energy correction problems which happens for soft or forward-backward jets. In 
addition initial and final state radiation effect is reduced by this cut. All jets have to be separated enough from the 
muon in the event by requiring 

Ai?(jet,,„uon) > 0.5 (2) 

Here Ai? = (Ayy)^ + (A^)^ and cj) is the azimuthal angle in the transverse plane. Figure |4] shows the number of 
reconstructed jets which passed the Et and \r]\ requirements mentioned above. Although this may not be a strongly 
distinguishing distribution, a cut on the number of jets is applied as the following: 

Number of jets > 3, with E^^?' > 30 GeV, < 2.5 (3) 

This cut ensures that events with low jet multiplicity (like single W events accompanied by 0, f , 2 or more jets) 
do not pass the signal selection cuts. As is seen from Fig. SI the number of reconstructed jets has a peak at three 
which is the result of reconstruction inefficiency and the cuts applied on the jet transverse energy and pseudorapidity. 
Appearance of events with more than four jets in Fig. 2] is a consequence of initial and final state radiation (ISR 
and FSR). They contribute to only ~ 14% of all events of the background and the signal with m{H^) = 80 GeV. 
The estimation of the error on the number of such events requires a reasonable understanding of the ISR and FSR 
performance. The additional jets may also be produced by experimental effects such as electronic noise or detector 
material activation and a detailed study of all such effects needs a full simulation of events at the presence of detector 
effects which is beyond the scope of this analysis. However as a reference, we mention an analysis of tt events in the 
lepton+jet final state which is similar to the final state studied in this work. This analysis, reported in (43| . uses 
different PYTHIA samples generated by varying parameters responsible for the amount of ISR and FSR radiation 
based on the approach described in [l0|. The conclusion is summarized in small uncertainties of about +1.7% 
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FIG. 4: Number of reconstructed jets passing Et > 30 GeV and |7;| < 2.5 in signal and background samples. 



and -1.3% on the measured tt cross section depending on increasing or decreasing the ISR and FSR effects. This 
uncertainty has been obtained for y/s — 7 TeV and should be evaluated for ^/s = 14 TeV, however, the overall 
conclusion is that we expect the uncertainty due to the extra jets associated with the tt pair to be at the level of few 
percent and it may not have a sizable effect on the signal significance. In the next step it is required that all jets 
should be separated by the AR requirement as the following: 



Ai?(jet.,jct,) > 0.5 (4) 

The indices i and j run over all selected jets in the event. Since at this point events contain three or more selected 
jets, this requirement makes sure that all pairs of those jets in the event are well separated and there is no overlap 
between the jet cones and they are well defined. 

The number of reconstructed and selected jets are reduced with heavier charged Higgs in the event as the energy of 
the b-jet produced in the top quark decay to charged Higgs is suppressed when a heavy charged Higgs is produced. 
Therefore the b-jet passes the kinematic cuts with a lower efficiency resulting in a lower jet multiplicity. This feature 
can also be verified by counting number of b-tagged jets which are identified if the reconstructed jet lies within a cone 
of AR < 0.3 around the parton level b-quark. An identification efficiency of 50% is also used when a jet matches a 
parton level b-quark. Figure [5] shows the number of b-tagged jets with the above algorithm. The cut on the number 
of b-tagged jets defined as: 



Number of b-tagged jets > 1, with £;^-taggcd jet ^ ^^^^ |,^b-taggcd jot| ^ 3.5 (5) 

not only selects tt events, but also rejects W-|-njets events which are produced with light jets, and also W'^W~ and 
ZZ events. The QCD multi-jet events are expected to be suppressed by this cut as well as the lepton requirement. 
In the next step, as there is a W boson in the event which decays to two light jets, a minimization procedure is 
performed to reconstruct the W boson and the top quark masses. The minimization is done by performing a loop 
over the jets in the event and finding the correct combination by minimizing the following x^- 



(TW C^top 

where the loop is over the jet indices I, m and n; mw and mtop are the nominal values of the W boson and the 
top quark masses, i.e. 80 and 175 GeV respectively and the ayy and atop values are taken to be 10 and 15 GeV 
respectively. Thanks to the kinematic cuts on jets this process is easily performed as every event has few hard jets to 
combine. Figures [Hand [7] show the reconstructed masses of W boson and top quark in signal events with m{H^) = 80 
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FIG. 5: Number of b-tagged jets in signal and background events. The b-tagging algorithm is based on matching jets which 
pass the kinematic cuts, with the parton level b-quarks. A b-tagging efficiency of 50% has been used. 
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FIG. 6: The reconstructed W candidate invariant mass in signal sample with m(_ff^) = 80 GeV. 

GeV. All other signal samples and the background sample show very similar distributions as the hadronic sector of 
these events, i.e. t — >■ W^h — >■ jjh is identical. 

To ensure that the selected event has a top quark and W boson candidate in the hadronic sector, a mass window 
cut is applied as the following: 

\mjj - mw\ < 10 GeV, \mjjj - mtop\ < 20 GeV (7) 

It should be noted that here no b-tagging has been used for the top reconstruction although in the previous step the 
b-jet requirement has already been applied. 

Finally the missing transverse energy is calculated as the negative sum over momentum components of all stable 
particles, which lie in the region \ri\ < 3.5, exclusing neutrinos in the event. A more general calculation should 
exclude sneutrinos, neutralinos, gravitons, gravitinos, however in this analysis the relevant sources of missing 
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FIG. 7: The reconstructed top quark candidate invariant mass in signal sample with m(H^) = 80 GeV. 
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FIG. 8: Missing transverse energy distributions in signal events with different charged Higgs massed and the background 
sample. 

Et are neutrinos. Of course in the real situation, this quantity is calculated from the imbalance of energy in the 
transverse plane of the detector. Figure |8] shows distributions of missing Et in signal and background. A cut on the 
missing Et in the event is applied as the following : 



Missing Et > 30 GeV 



(8) 



As is seen from Fig. |8l distributions are similar producing no sizable increase in the signal significance when harder 
cuts are applied. 

Another quantity which can be made is the charged Higgs candidate transverse mass. This is calculated as the 
following: 



TUT = ^2.pr°"-^^T'"-(l - C0SA<^(„,„„,^™)) 



(9) 
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FIG. 9: The charged Higgs candidate transverse mass distribution in signal samples with different charged Higgs massed and 
the background sample. 





80 GeV 


100 GeV 


120 GeV 


150 GeV 


o-x BR 


302 fb 


228 fb 


149 fb 


38 fb 


Muon selection (Eq. [T]) 


72.1 % 


79.6 % 


85.5 % 


91% 


Number of jets > 3 (Eq. E]) 


72.4 % 


70.2 % 


65.9 % 


52.5 % 


AR{j,j) > 0.5 (Eq.H 


87.6% 


87.2% 


87% 


88% 


Number of b-jets > 1 (Eq. [5]) 


62.7 % 


61.9 % 


59.7 % 


49.7% 


W and top mass window (Eq. [T]) 


79.1 % 


79.2 % 


79.9% 


80.6% 


MET > 30 GeV (Eq. [8]) 


69.1 % 


78.4 % 


85% 


91.2% 


Charged Higgs mass window (Eq. 1101) 


70.7 % 


71.8% 


74.4% 


78.9 % 


Total efficiency 


11.1 % 


13.47% 


14.8 % 


12.1% 


Events @ 300 fb'^ 


10057 


9213 


6616 


1379 



TABLE II: Signal selection efficiencies. 



where Acji/^^^^^^ £;™'==) is the angle between the muon flight and the missing Et in the transverse plane. Figure [3] 
shows the distribution of this quantity in signal samples and the background. The lower cut applied on the charged 
Higgs transverse mass is set as the following: 



Charged Higgs candidate tot > 50 GeV (10) 

As the signal to background ratio is small, harder cuts end up with smaller statistics while not increasing the signal 
significance sizably. Now the signal and background relative and total efficiencies are calculated to evaluate the signal 
significance when all cuts are applied. Table HIl shows the signal selection efficiencies for different charged Higgs masses. 
Table Hill shows the background efficiencies and the final event yield. 

The numbers quoted in Tabs. HI] and IHII lead to the following results listed in Tab. IIVI As is seen the signal 
significance is close to 5cr near the direct search lower limit. Although the signal efficiency increase with higher 
charged Higgs masses which is due to the harder kinematics of events, the signal significance decreases due to the 
reduction in cross section. A problem to note is that the signal to background ratio is very small making the statistical 
conclusions about the existence of the signal, a challenging task. 
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4315000 



TABLE III: tt background selection efficiencies. 





80 GeV 


100 GeV 


120 GeV 


150 GeV 


S/B 


0.002 


0.002 


0.001 


0.0003 


s/Vb 


4.8 


4.4 


3.2 


0.66 



TABLE IV: Signal to Background ratio and the signal significance. 

IV. DISCUSSION ON THE ACCURACY AND CONSISTENCY OF RESULTS 

The analysis presented here does not reflect completely the LHC experimental environment. A detailed study of 
this channel requires a full simulation of an LHC detector including all experimental environment in detail followed 
by a full analysis of experimental and theoretical uncertainties involved in the analysis which is beyond the scope 
of this paper. This can be done by one of LHC experiments in future with access to the full set of a real analysis 
software. In this work, however, a result consistency check is performed between different software versions. 
The PYTHIA 6.4.21 used for the analysis is one of the most recent versions of the 6*'' series of this package. Although 
PYTHIA 8 series are taking over, PYTHIA 6 is still being used by LHC experiments in their physics simulations. 
Therefore we rely on PYTHIA 6 results as the basis. 

Concerning HDECAY, the up-to-date version 4.1 was compared with the version used in the analysis, i.e., version 3.5. 
No sizable difference between the branching ratios given by the two above versions was observed. 
The MCFM version 5.7 which was used in the analysis was also compared with the up-to-date version 6.1. In this 
case also, no change in the cross sections was observed up to the femtobarn level. 

The choice of the PDF set can potentially have some impact on the event cross section as well as selection efficiencies. 
This was verified by using LHAPDF 5.8.6 (the up-to-date version) and using three sets of PDFs, i.e., MSTW (2008), 
CTEQ6mE (2002) and NNPDF 2.1 with 100 fit parameters (2011). The MSTW and MRST package used in the 
analysis calculate as evolution up to NNLO in its perturbation series, while CTEQ6mE and NNPDF 2.1 evolve as 
up to NLO. Results for signal and backgrounds are listed in Tabs. IVland lvTl As seen from these tables, the overall 
change in the final results in terms of signal significance is expected to be at the level of few percent. There can 
be more effects due to as variations in PDFs which was not studied in this analysis, however, the overall effect of 
PDF-t-Qfs uncertainty in the results is expected to remain at the level of few percent. An example of such analyses, 
although done for a different process, can be found in 44]. 

Finally, consistency of the results can be verified through Figs. I3l4l5l8l and[^ which reveal a general effect of increasing 
the charged Higgs mass which produces harder events with higher muon px, missing transverse energy, etc. On the 
contrary the heavier the charged Higgs, the less number of jets identified. This is due to the fact that the 6-jet 
produced in the top quark decay is suppressed when increasing the charged Higgs mass resulting in a softer jet which 
is identified with a smaller probability as it may not pass the transverse energy threshold applied in the analysis. 

V. CONCLUSIONS 

The charged Higgs decay to muon and muon neutrino was studied in this work. Results show that the light charged 
Higgs may in fact be observable through this decay channel. The signal significance ranges between 4.8ct to 0.66(t 
for 80 GeV < m{H^) < 150 GeV with tan/3 = 20. The signal yield is few thousand events, however, a challenging 
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PDF set 


efHciency 


Cross section 


MSTW 2008 


-2% 


0.2% 


CTEQ6mE 


-0.3% 


-2% 


NNPDF2.1 


4% 


2% 



TABLE V: Relative change in the total efficiency and cross section of the signal compared to the PDF set used in the analysis 
(MRST 2006) 



PDF set 


efficiency 


Cross section 


MSTW 2008 


4% 


-0.1% 


CTEQ6mE 


-1% 


-1% 


NNPDF2.1 


1% 


4% 



TABLE VL Relative change in the total efficiency and cross section of the background compared to the PDF set used in the 
analysis (MRST 2006) 

problem would be the very Small S/B ratio which is of the order of 10~^ or smaller. The singal significance decreases 

rapidly with increasing the charged Higgs mass, however, the higher charged Higgs mass, the smaller the excluded 
area in {m{H^), tan/3 ) plane. As an example, for m{H^) > 150 GeV, tan/3 values are not excluded up to tan/3 ~ 50. 
This point in the parameter space leads to a signal significance of 3.7a. Although other charged Higgs decay channels 
produce higher signal rates, their 5(7 contours do not cover the parameter space completely. Therefore adding this 
channel with a proper statistical method may increase the overall sensitivity of LHC detectors to the charged Higgs 
signal. 
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